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ABSTRACT 

It  is  necessary  to  clarify  the  effect  of  Marangoni  convection  on  single  crystal  growth  from  a 
melt  in  order  to  improve  the  quality  of  the  grown  crystal.  Particularly,  the  deviation  of  crystal- 
melt  (C-M)  interface  from  a planar  shape  is  a major  problem  because  it  may  deteriorate  the  quality 
of  the  grown  crystal.  In  this  paper,  we  investigated  the  effect  of  thermal  and  solutal  Marangoni 
convection  on  C-M  interface  shape  in  an  In-Sb  binary  system  by  the  horizontal  Bridgman  (HB) 
method.  The  C-M  interface  concavity  strongly  depends  on  the  cooling  rate  and  the  temperature 
gradient  under  uniform  concentration  distribution  conditions  in  the  melt.  A large  concavity  was 
observed  at  low  cooling  rates  and  high  temperature  gradient  conditions.  The  concavity  was  found 
to  be  caused  by  thermal  Marangoni  convection,  by  taking  Peclet  number  into  account.  Then,  we 
varied  the  composition  of  the  In-Sb  binary  system  to  induce  solutal  Marangoni  convection 
intentionally.  The  C-M  interface  was  kept  planar  in  case  solutal  Marangoni  convection  occurred  in 
the  direction  opposite  to  the  thermal  one.  Therefore,  we  believe  that  the  utilization  of  solutal 
Marangoni  convection  will  be  a new  control  technique  to  make  the  C-M  interface  planar  for  the 
HB  system.  From  these  results,  it  was  clarified  that  Marangoni  convection  plays  a significant  role 
in  the  HB  crystal  growth  system, 

INTRODUCTION 

At  present,  most  of  the  semiconductor  bulk  crystals  are  grown  by  melt  growth  methods  such 
as  the  Czochralski  (CZ)  or  the  horizontal  Bridgman  (HB)  method.  During  crystal  growth  by  these 
methods,  convection  in  melt  affects  the  quality  of  the  grown  crystals.  Therefore,  it  is  necessary  to 
clarify  the  mechanism  of  convection  in  melt  and  to  control  the  convection  to  improve  the  quality 
of  the  crystal.  The  convection  can  be  divided  into  buoyancy  convection  and  Marangoni 
convection.  The  driving  force  of  the  former  is  the  buoyancy  difference  in  the  melt,  and  that  of  the 
latter  is  the  surface  tension  difference  on  the  free  surface.  The  convection  can  also  be  divided  into 
thermal  and  solutal  convection  because  those  driving  forces  vary  with  temperature  and 
concentration.  Particularly,  thermal  buoyancy  convection  has  been  extensively  investigated  since 
as  early  as  the  beginning  of  the  20th  century  [1].  On  the  other  hand,  the  importance  of  Marangoni 
convection,  which  has  nothing  to  do  with  gravity,  was  recognized  after  semiconductor  crystal 
growth  experiments  were  conducted  under  microgravity  conditions  in  the  1980s  [2].  Since  then, 
Marangoni  convection  has  been  investigated  by  many  researchers  [3,4]. 

Regarding  solutal  convection,  although  solutal  buoyancy  convection  has  been  studied  as  a 
double  diffusive  problem,  solutal  Marangoni  convection  had  not  been  studied  so  far.  Therefore, 
we  performed  some  experiments  to  measure  surface  velocities  caused  by  buoyancy  convection 
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and  Marangoni  convection  using  a rectangular  boat  with  a free  surface.  Consequently,  we  found 
that  the  surface  velocities  of  solutal  Marangoni  convection  arc  approximately  3-5  times  higher 
than  those  of  thermal  Marangoni  convection  [5],  Moreover,  in  the  ca.se  of  coexistence  of  thermal 
and  solutal  Marangoni  convection,  solutal  Marangoni  convection  was  found  to  be  dominant  in  the 
surface  flow,  and  the  flow  direction  could  be  opposite  to  that  without  solutal  convection  [6j. 

Also,  we  carried  out  parabolic  flight  experiments  and  found  that  Marangoni  convection  is 
dominant  in  the  surface  flow  in  case  both  buoyancy  convection  and  Marangoni  convection  coexist 
[6].  From  these  re.sults,  we  demonstrated  that  Marangoni  convection,  especially  solutal 
Marangoni  convection,  is  significant  in  a rectangular  boat  such  as  that  u.sed  in  the  HB  method.  As 
one  effect  of  Marangoni  convection  on  the  grown  crystals,  the  change  of  C-M  interface  shape  is 
often  discussed.  The  deviation  of  the  C-M  interface  shape  from  a planar  one  is  a major  problem 
because  it  is  assumed  to  lower  the  quality  of  the  grown  crystals.  Using  the  HB  method,  Lan  et  al. 
clarified  numerically  [7]  and  experimentally  [8]  that  the  effect  of  Marangoni  convection  on  C-M 
interface  shape  is  significant. 

In  this  study,  we  performed  InSb  crystal  growth  experiments  by  the  HB  method  to 
investigate  the  effects  of  thermal  and  .solutal  Marangoni  convection  on  C-M  interface  shape.  We 
varied  the  composition  of  the  In-Sb  binary  sy.stem  to  induce  .solutal  Marangoni  convection 
intentionally  in  addition  to  thermal  Marangoni  convection.  Ba.scd  on  the  results,  we  sugge.sted  a 
simple  way  of  controlling  C-M  interface  shape. 

EXPERIMENT 
Concept  of  experiment 

First,  we  discuss  the  interaction  between  thermal  and  solutal  Marangoni  convection.  As 
described  in  the  Introduction,  Marangoni  convection  is  induced  by  surface  tension  difference  due 
to  the  difference  of  temperature  or  concentration  on  a free  surface.  In  our  previous  study  [6],  it 
was  clarified  that  the  flow  pattern  of  Marangoni  convection  could  be  divided  into  three  cases  as 
follows.  Here,  x denotes  the  composition  ofSb  in  the  In-Sb  binary  system  (lni.,Sb,). 

(a)  Normal  ca.se  (x=0.5) 

In  this  case,  there  is  no  concentration  difference  in  the  melt.  Flow  of  thermal  Marangoni 
convection  is  induced  in  the  direction  shown  in  figure  1(a).  This  is  because  surface  tension  a near 
the  C-M  interface  is  larger  than  that  at  the  bulk  of  the  melt. 

(b)  Acceleration  case  (0<x<0.5) 

Thermal  Marangoni  convection  is  induced  in  the  same  direction  as  (a).When  crystal  growth 
proceeds,  concentration  of  Indium  Ci„  near  the  C-M  interface  becomes  higher  than  that  at  the  bulk 
of  the  melt.  Because  it  is  known  that  Oi„  is  larger  than  Ost,,  surface  tension  near  the  C-M  interface 
is  larger  than  that  at  the  bulk  of  the  melt.  This  causes  solutal  Marangoni  convection.  Therefore, 
the  flow  direction  of  solutal  convection  is  the  same  as  that  of  thermal  Marangoni  convection. 

(c)  Deceleration  ca.se  (0.5<x<0.68) 

On  the  other  hand,  because  C|„  near  the  C-M  interface  becomes  lower  than  that  at  the  bulk  of 
the  melt,  surface  tension  at  the  bulk  of  the  melt  becomes  larger  than  that  near  the  C-M  interface. 
Therefore,  the  flow  of  .solutal  Marangoni  convection  is  induced  in  the  direction  opposite  to  that  of 
thermal  Marangoni  convection.  Because  it  was  clear  that  solutal  Marangoni  convection  is 
stronger  than  thermal  Marangoni  convection  [5],  the  flow  direction  becomes  opposite  to  that 
without  considering  solutal  Marangoni  convection. 
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(a)  Normal  case 
(x=0.5) 


In  rich 


cold  ^ hot 


(b)  Acceleration  case 
(0<x<0.5) 


Sb  rich 


cold  ^ ^ hot 

(c)  Deceleration  case 
(0.5<x<0.68) 


; Thermal  Marangoni  convection 
: Solutal  Marangoni  convection 


Figure  1.  Flow  pattern  by  interaction  between  thermal  and  solutal  Marangoni  convection 


Experimental  procedure 


To  avoid  contamination  and  formation  of  metal  oxide  film  on  the  free  surface  of  the  melt, 
special  precautions  were  taken.  In  preparing  the  test  sample,  Ini.xSb,,  In  and  Sb  (5N, 
manufactured  by  Dowa  Mining)  were  mixed  at  the  desired  ratio,  and  the  mixture  was  dissolved 
and  quenched  in  an  oxygen  reducing  atmosphere  (Ar97%+H23%).  Then,  chemical  etching  of  the 
test  sample  was  carried  out.  Moreover,  a carbon  crucible  was  baked  for  degassing.  Figure  2 
shows  the  cross  section  of  the  carbon  crucible.  The  test  sample  was  set  in  the  reservoir  tank,  and 
InSb,  which  has  a length  of  approximately  10mm,  was  also  set  at  the  cold  side  of  the  observation 
part  as  a seed  crystal.  After  heating  the  crucible  and  the  reservoir  tank,  the  melted  sample  flowed 
into  the  observation  part  through  a small  hole  and  contacted  with  the  seed  crystal.  The  small  hole 
prevents  efflux  of  the  remaining  metal  oxides  in  the  sample.  The  experimental  system  is  shown  in 
figure  3.  The  temperatures  of  the  hot  side,  Th,  the  cold  side,  Tc,  and  the  middle  of  the  crucible, 

Tm  were  controlled  by  PID  temperature  controllers  and  power  controllers.  The  temperature 
difference  between  Th  and  Tc,  and  cooling  rates  of  Tn,  Tm,  and  Tc  were  kept  constant  as  shown 
in  table  1 . To  reduce  buoyancy  convection,  we  used  a shallow  crucible.  The  length,  width  and 
depth  of  the  cavity  of  the  crucible  were  48,  10,  and  5 mm,  respectively.  The  test  sample  was 
solidified  unidirectionally  in  the  crucible.  C-M  interface  shape  was  observed  by  means  of  a 3-CCD 
video  camera.  C-M  interface  shape  was  evaluated  in  terms  of  the  degree  of  concavity  (figure  4). 
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Figure  4.  Top  view  of  C-M  interface  shape 
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RESULTS  AND  DISCUSSION 


We  varied  the  control  parameters  such  as 
cooling  rate,  temperature  gradient,  and 
composition  of  In-Sb  in  the  experiments.  The 
experimental  conditions  are  shown  in  table  1 . 

Figure  5 shows  concavity  as  a function  of 
Lf/Lo  at  various  cooling  rates  with  the  same 
temperature  gradient  of  14.rC/cm.  Here,  Lc 
and  L<)  denote  the  length  of  the  crystal  during 
growth,  and  the  total  crystal  length  after 
growth,  respectively.  Concavity  increa.sed  as 
the  growth  proceeded,  and  then  saturated  at 
an  Lf/L(i  of  0.5.  Moreover,  the  C-M  interface 
concavity  obtained  after  the  growth  at  low  cooling  rates  was  larger  than  that  at  high  cooling  rates. 
This  is  because  long-term  exposure  to  thermal  Marangoni  convection  at  low  cooling  rates  resulted 
in  the  increase  of  the  erosion  of  the  C-M  interface.  Figure  6 shows  concavity  as  a function  of 
temperature  gradient.  The  temperature  gradient  was  varied  from  8.4  to  I4.1°C  /cm  at  the  .same 
cooling  rate  of-5°C/h.  The  concavity  was  found  to  depend  strongly  on  the  temperature  gradient. 
Although  in  figure  5.  the  concavity  under  the  temperature  gradient  of  -S^C/h  was  approximately 
0.55,  the  concavity  was  kept  0 at  8.4°C/cm  even  if  the  cooling  rate  w;is  -5®C/h  in  figure  6.  The 
low  temperature  gradient  indicates  that  thermal  Marangoni  convection  is  weak.  From  figures  5 
and  6,  there  seems  to  be  a strong  correlation  between  driving  force  of  Marangoni  convection  and 
C-M  interface  shape. 


Table  I.  Experimental  conditions 


Run  No. 

Cooling  rate 

rc/hi 

Temperature  gradient 

Composition  x 
of  Inl-xSbx 

1 

-5 

14.1 

0.5 

2 

-M.7 

14.1 

0.5 

3 

-16.7 

14.1 

0.5 

4 

-19.7 

14.1 

0.5 

5 

-5 

11.9 

0.5 

6 

-5 

8.4 

0.5 

7 

-16.7  i 

14.1 

0.4 

8 

-16.7 

14.1 

0.6 

Figure  5.  Concavity  as  a function  of  Lc/Ui 
The  dashed  line  indicates  the  saturation 
point  of  an  increase  of  concavity. 


Figure  6.  Concavity  as  a function  of 
temperature  gradient 


In  order  to  confirm  this  correlation,  we  determined  Marangoni  number  (Ma)  from 
temperature  profiles  and  C-M  interface  position.  Ma  represents  the  driving  force  of  Marangoni 
convection.  Ma  is  defined  by  eq.(  1 ). 


Ma  = 


HV 


eq.(l) 
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In  our  previous  study  [9],  we  found  a relationship  between  Ma,  Pr,  and  Re,  as  shown  in 

eQ.(2). 

Re  = 0.28  ■ (Ma  eq.(2) 

Here,  Re  (Reynolds  number)  and  Pr  (Prandtl  number)  represent  flow  state  of  fluids  and  heat 
transport  phenomena,  respectively.  Pr  is  defined  by  eq.(3). 

Pr=-  eq.(3) 

a 

From  eqs.(l),  (2),  and  (3),  we  were  able  to  estimate  Re.  After  estimation  of  Re,  we 
determined  Peclet  number  (Pe)  using  eq.(4). 

Pe  = Pr-  Re  eq.(4) 

Pe  indicates  whether  convection  or  conduction  dominates  the  heat  transport  phenomena. 
Therefore,  Pe  indicates  the  effect  of  Marangoni  convection  on  the  heat  transfer  through  the  C-M 
interface. 

Figure  7 shows  the  calculated  Pe  as  a function  of  Lc/Lo.  Pe  decreased  gradually  as  growth 
proceeded  up  to  an  Lc/Lo  of  0.5,  and  then  abruptly  decreased  thereafter.  In  other  words,  the 
effect  of  thermal  Marangoni  convection  on  the  heat  transfer  decreased  abruptly  at  Lc/Lo  =0.5. 
This  value  agrees  with  the  saturation  point  of  the  increase  of  concavity  shown  in  figure  5.  From 
these  results,  it  was  clarified  that  concavity  is  mainly  due  to  Marangoni  convection. 

As  a simple  way  to  keep  C-M  interface  planar,  we  suggest  that  the  Marangoni  convection  be 
restrained.  The  flow  direction  may  be  opposite  to  that  of  the  normal  case  in  the  presence  of 
solutal  Marangoni  convection,  as  reported  in  our  previous  study  [6].  Therefore,  we  varied  the 
composition  of  the  In-Sb  binary  system  to  induce  solutal  Marangoni  convection  intentionally  in 
addition  to  thermal  Marangoni  convection.  Figure  8 shows  the  concavity  as  a function  of  L</Lo. 
Although  the  increase  of  concavity  in  the  case  of  x=0.4  was  larger  than  that  in  the  case  of  x=0.5, 
increase  in  the  case  of  x=0.6  was  much  .smaller  than  that  in  the  case  of  x=0.5.  This  is  becau.se  flow 
velocity  was  accelerated  in  the  case  of  x=0.4  and  decelerated  in  the  case  of  x=0.6.  Moreover,  we 
could  keep  the  C-M  interface  planar  in  the  case  of  x=0.6.  From  these  results,  we  found  that  we 
can  keep  a planar  interface  by  inducing  solutal  Marangoni  convection  intentionally  in  the  direction 
opposite  to  that  of  thermal  Marangoni  convection. 


Figure  7.  Pe  as  a function  of  Lc/Lo  Figure  8.  Concavity  as  a function  of  Lc/Lo 
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CONCLUSIONS 


We  investigated  the  effect  of  thermal  and  solutal  Marangoni  convection  on  crystal-melt  (C- 
M)  interface  shape  in  the  In-Sb  binary  system  by  the  horizontal  Bridgman  (HB)  method.  The  C-M 
interface  concavity  strongly  depends  on  the  cooling  rate  and  the  temperature  gradient  under 
uniform  concentration  distribution  conditions  in  the  melt.  A large  concavity  was  observed  at  low 
cooling  rates  and  high  temperature  gradient  conditions.  The  concavity  was  found  to  be  due  to 
thermal  Marangoni  convection,  taking  Peclet  number  into  account.  We  varied  the  composition  of 
the  In-Sb  binary  .system  to  induce  solutal  Marangoni  convection  intentionally.  The  C-M  interface 
remained  planar  when  solutal  Marangoni  convection  occurred  in  the  direction  opposite  to  the 
thermal  one.  Therefore,  we  believe  that  the  utilization  of  solutal  Marangoni  convection  will  be  a 
new  technique  to  control  the  planarity  of  the  C-M  interface  in  the  HB  system.  From  these  results, 


it  was  clarified  that  Marangoni  convection  plays 
system. 

NOMENCLATURE 

Cin  : concentration  of  Indium  [at%] 

L)  ; total  length  of  crystal  [m] 

Lc  : length  of  crystal  during  growth  [m] 
Lii  : length  of  melt  [m] 

Ma  : Marangoni  number  [-] 

Pe  : Peclet  number  [-] 

Pr  ; Prandtl  number  [-] 

Re  : Reynolds  number  [-] 

T : temperature  [°C] 

Tr  : temperature  at  the  cold  side  [°C] 

Tii  : temperature  at  the  hot  side  [®C] 

Tm  : temperature  at  the  middle  of  the 

crucible  [°C] 

h ; deviation  of  crystal-melt  interface  [m 

d ; width  of  melt  [m] 
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AT  ; temperature  difference  [K] 

a : heat  diffusion  coefficient  [m"s  '] 

V : kinematic  viscosity  [mV] 

p : viscosity  [kguV's"'] 

p ; density  [kgnf’] 

a ; surface  tension  [kg.s"'] 

ai„  ; surface  tension  of  Indium  [kgs'"] 

asb  : surface  tension  of  Antimony  [kgs‘‘] 
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